Bifurcations and chaos in semiconductor superlattices with a tilted magnetic field. 
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We study the effects of dissipation on electron transport in a semiconductor superlattice with an 
applied bias voltage and a magnetic field that is tilted relative to the superlattice axis. In previous 
work, we showed that although the applied fields are stationary, they act like a THz plane wave, 
which strongly couples the Bloch and cyclotron motion of electrons within the lowest miniband. As 
a consequence, the electrons exhibit a unique type of Hamiltonian chaos, which creates an intricate 
mesh of conduction channels (a stochastic web) in phase space, leading to a large resonant increase 
in the current flow at critical values of the applied voltage. This phase-space patterning provides a 
sensitive mechanism for controlling electrical resistance. In this paper, we investigate the effects of 
dissipation on the electron dynamics by modifying the semiclassical equations of motion to include a 
linear damping term. We demonstrate that even in the presence of dissipation, deterministic chaos 
plays an important role in the electron transport process. We identify mechanisms for the onset 
of chaos and explore the associated sequence of bifurcations in the electron trajectories. When the 
Bloch and cyclotron frequencies are commensurate, complex multistability phenomena occur in the 
system. In particular, for fixed values of the control parameters several distinct stable regimes can 
coexist, each corresponding to different initial conditions. We show that this multistability has clear, 
experimentally-observable, signatures in the electron transport characteristics. 

PACS numbers: 05.45.-a, 72.20.Ht, 73.21.Cd 



I. INTRODUCTION. 



Semiconductor superlattices (SLs) are nanostructures 
made from alternating layers of two different semicon- 
ductor materials, usually with very similar lattice con- 
stants, for example GaAs and (AlGa)As Due to the 
different energy band gaps of the two materials, the con- 
duction band edge of an ideal superlattice is periodically 
modulated. Typically, a SL made by molecular beam 
epitaxy contains 10-100 quantum wells in series, which 
are coupled by tunnel barriers. This periodic potential 
leads to the formation of energy bands, known as "mini- 
bands" , for electron motion perpendicular to the layers 
A voltage applied to ohmic contacts at the two ends 
of the SL generates an electric field, F, perpendicular to 
the plane of the layers (see Fig. [1]), which causes charge 
to flow through the device. The current-voltage char- 
acteristics of SLs are usually highly nonlinear due to a 
variety of quantum mechanical effects including resonant 
tunneling, the formation of Wannier-Stark energy level 
ladders, and the occurrence of Bloch oscillations [l], 
whose frequency is proportional to the spatial period, 
d, of the SL and, also, to F. In natural crystals, d is 
so small (~ 0.3 nm) that Bloch oscillations do not oc- 
cur because the Bloch frequency is much less than the 
electron scattering rate. But in SLs, d and the corre- 
sponding Bloch frequency can be large enough for the 
electrons to perform Bloch oscillations, which then play 
a key role in both the dc and high-frequency charge trans- 
port processes. The onset of Bloch oscillations localizes 
the electrons, thus causing their drift (average) velocity 
to decrease as the electric field increases. This negative 
differential velocity can induce high-frequency collective 
oscillations of the conduction electrons within the SL lay- 
ers [3] , making SLs attractive for the generation and de- 




FIG. 1: Schematic layer structure of a semiconductor SL 
formed from two different semiconductor materials, shaded 
light and dark gray. The co-ordinate axes show the orien- 
tation of the tilted magnetic field, B, which lies in the x — z 
plane at an angle 6 to the SL (x) axis. The lower arrow shows 
direction of the electric field, F, applied perpendicular to the 
plane of the layers and anti-parallel to the x axis. 



tection of electromagnetic radiation in the GHz to THz 
frequency range [J, 0, IS ■ 

The occurrence of negative differential velocity for elec- 
trons in SLs has also led to fundamentally new regimes 
of charge transport involving deterministic chaos, charac- 
terized by complex irregular electron dynamics [2, H, 11 • 
Understanding the effects of chaos on the current-voltage 
characteristics, I(V), and high-frequency electromag- 
netic properties of SLs is an emerging research area at 
the interface between nonlinear dynamics and condensed 



2 



matter physics. In most previous studies, chaos in SLs 
was identified in the cooperative motion of interacting 
electrons in either periodically driven [HI, HH or undriven 
SLs Q. 

Recent theoretical and experimental work has revealed 
that the single-particle trajectories and collective behav- 
ior of electrons moving through the lowest miniband of 
a biased SL with a tilted magnetic field have unique 
properties [HI, EH, EE EBj- For example, the electrons 
exhibit an unusual type of Hamiltonian chaos, known 
as "non-KAM" chaos (l7l [H, [l9[ , which does not obey 
the Kolmogorov-Arnold-Moser (KAM) theorem and pro- 
vides a sensitive new mechanism for controlling electron 
transport. Remarkably, the effective classical Hamilto- 
nian for the electron motion in a SL has an intrinsi- 
cally quantum-mechanical origin as it depends explicitly 
on the energy versus wavevector dispersion relation of 
the miniband [H, [3] ■ Analysis of Hamilton's equations 
reveals that the stochastic electron motion switches on 
abruptly when the field parameters satisfy certain reso- 
nance conditions, described in detail below. The onset 
of chaos delocalizes the electrons by imprinting an intri- 
cate mesh of conduction channels (known as a "stochas- 
tic web") in phase space, which produces large resonant 
enhancement of the electron velocity and current flow 
measured in experiment This phase-space pattern- 
ing provides a fundamentally new concept for controlling 
electrical conductance [2(| and operates even at room 
temperature 

In the absence of dissipation, the resonant derealiza- 
tion of the conduction electrons would produce a series 
of delta-function peaks in the /(V) curves. But in real 
SLs, the resonances are broadened because the electrons 
scatter both elastically [ij, due mainly to roughness at 
the interfaces between the tunnel barriers and quantum 
wells and to ionized donor atoms, and inelastically via 
the emission and absorption of phonons. When calculat- 
ing the transport characteristics of the SLs, the effect of 
scattering on the electron dynamics must be taken into 
account in order to obtain I(V) curves that agree with 
experiment. In our previous theoretical work (T3 , fl4| . 
we incorporated electron scattering in the following phe- 
nomenological way. First, we calculated the electron tra- 
jectories by solving Hamilton's equations for the system, 
assuming no scattering. Next, we used these trajectories 
to determine the average electron velocity from a simple 
kinetic formula in which scattering appears a posteriori 
as an exponential damping term whose decay rate de- 
pends on the elastic and inelastic relaxation times. In 
this paper, we investigate whether the striking electron 
resonance effects and phase space patterning predicted by 
our earlier collision-free Hamiltonian model of the elec- 
tron dynamics persist when scattering is included a priori 
in the equations of motion. We find that when the scat- 
tering is described in this way, it creates far richer elec- 
tron dynamics than expected from our previous Hamil- 
tonian analysis. In particular, scattering makes the sys- 
tem highly sensitive to changes in the control parameters 



(electric and magnetic fields), which induce complex bi- 
furcation sequences characterized by alternating windows 
of stability and chaos. Counter-intuitively, dissipation 
also enhances the resonant derealization of the semi- 
classical trajectories by creating attractors that drive the 
electrons rapidly through the SL. The regimes of dissi- 
pative chaos that we have identified have a pronounced 
effect on the velocity- field characteristics of the mini- 
band electrons and provide mechanisms for controlling 
SL magnetotransport by exploiting complex nonlinear 
dynamics. 

The structure of the paper is as follows. In Section 
II, we introduce the semiclassical equations of motion for 
a miniband electron in a tilted magnetic field and show 
how scattering is included a priori in those equations. In 
Section III, we show that scattering dramatically changes 
the electron trajectories and phase space structure of the 
system by creating attractors corresponding to periodic 
orbits that extend through the whole SL. In Section IV, 
we present a detailed analysis of the stability of the elec- 
tron orbits. Section V explores the complex series of res- 
onances that scattering induces in the electron dynam- 
ics and elucidates the effect of these resonances on the 
electron velocity and experimentally-measured electrical 
current. Finally, in Section VI, we summarize our results 
and draw conclusions. 



II. MODEL EQUATIONS. 

In a SL, quantum mechanical tunneling broadens the 
discrete quantized energy levels of each individual quan- 
tum well into SL minibands. The miniband states are 
delocalized Bloch functions specified by the crystal mo- 
mentum p = (p x ,p y ,p z )=fik, where k is the correspond- 
ing electron wavevector. Within the tight-binding ap- 
proximation, the energy versus crystal momentum dis- 
persion relation for the lowest miniband is E(p) = A[l — 
cos(p x d/H)]/2 + (py+p1)/2m*, where A is the miniband 
width, d is the SL period, and m* is the electron ef- 
fective mass for motion in the y — z plane. The crystal 
momentum component, p x , is taken to lie within the first 
minizone of the SL. 

Throughout this paper, we consider electron motion in 
an electric field F = (— F , 0, 0) applied anti-parallel to the 
x axis, and a tilted magnetic field B = (B cos 9,0, B sin 6) 
(Fig. [T]). In a semiclassical picture, which neglects inter- 
miniband tunneling, the force produced by the electric 
and magnetic fields changes the electron's crystal mo- 
mentum at a rate 

^ = - e (F+(V pJ B(p)xB), (1) 

where e is the electronic charge. Equation (I) can be 
written in the component form 

p x = eF — uj c py tan 9 (2) 
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UcPz 



(3) 
(4) 



where the left hand terms are time derivatives of the 
crystal momentum components and Co c —eB cos 9 /m* is 
the cyclotron frequency corresponding to the magnetic 
field component along the x axis. It follows from Eqs. 
(2-4) that 



Pz 



C s\a{Kp z — uigt + • 



(5) 



where C = (-m*u 2 dAtaii9) /2h, K = dt&n9/H, and 
u>b = eFd/H is the Bloch frequency. The phase, 4> — 
d(p x (t = 0) + p z (t — 0)taxi9)/H, depends on the initial 
conditions and equals zero for electrons starting from rest 
EH]- Equation (5) describes completely the electron 
motion because its solution, p z (t), uniquely determines 
all of the other dynamical variables [l4| ■ Consequently, to 
include the effects of dissipation ab initio in the equations 
of motion, we formally introduce a relaxation term, ap z , 
into the left hand side of Eq. (5), which becomes 



Pz 



ap z + to 2 c p z = C sin(Kp z - uj B t 



(6) 



where a is the damping constant. In this approach, a 
has the same meaning as a constant collision rate in the 
Boltzmann equation. The dissipation term in Eq. ^ 
directly affects the in-plane momentum component, p z . 
But when 9^0, electron motion in the x, y, and z direc- 
tions is coupled, which means that the dissipation term 
affects all of the dynamical variables, thus simulating the 
scattering processes that occur in a real SL device. 

Recent experiments [HI, [22| have shown that the pres- 
ence of a high magnetic field strongly affects the inelas- 
tic scattering mechanisms required for current to flow 
through a biased SL [23J. In particular, when 9 = 
and huj c > A, the emission of longitudinal optical (LO) 
phonons is strongly suppressed if the phonon energy, 
fkdLo — 36 meV in GaAs, exceeds A [U [22|, which is 
the case for the SL that we studied in Rcf. [l4j]. In 
this low-dissipation regime, the current falls almost to 
zero because the oscillatory Bloch motion along the x 
direction is not damped by LO phonon emission. Our 
model equations (2-6) capture this magneto-suppression 
of the current, because Eqs. (2), (4), and (6) decou- 
ple in the limit 9 — > 0, which means that the x motion 
is undamped, and so the mean electron velocity is zero. 
In the experiments [l4j, as 9 increases from 0, the cur- 
rent rises dramatically because the selection rules that 
prevent LO phonon emission when 9 = are broken. 
Our model equations (2-6) describe this activation of LO 
phonon emission by increasing the coupling between mo- 
tion parallel and perpendicular to x as 9 increases from 0, 
thus enhancing scattering-assisted transport through the 
SL. Consequently, the equations encapsulate the strongly 
^-dependent scattering rates observed experimentally in 
the high magnetic field regime that is the focus of this pa- 
per. As a consequence, however, they do not accurately 



describe the strongly damped Bloch motion that occurs 
in the limit B = 0, which we do not consider here. 

We used Eq. (|6|) to investigate the effect of dissipa- 
tion on the motion of electrons in SLs similar to those 
described and studied experimentally in [l3l . Il4j , taking 
A = 26.2 meV, d = 10.3 nm, and m* = 0.067m e , where 
m e is the mass of a free electron. In our analysis, we 
first solve Eq. (6) numerically to obtain p z (t), which we 
then use to determine the other crystal momentum com- 
ponents 



Px 



Pz 

Py = r~, 



p x {t = 0) + eFt - (p z - Pz (t = 0)) tan 9, 



and the electron velocity components 



x = — sm(Kp z - uj B t 



V = 



Pz 



ui c mr 



(7) 
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Eq. ([6]) describes the motion of a damped harmonic 
oscillator, whose natural frequency equals the cyclotron 
frequency u> c , driven by a plane wave of wave number 
K and frequency equal to the Bloch frequency u>b- The 
linear dissipation term describes the coupling of the oscil- 
lator to a "reservoir" of many additional degrees of free- 
dom, which, for the particular case of electrons in a SL, 
is provided by elastic and inelastic scattering processes. 

An undamped harmonic oscillator driven by a plane 
wave is one of few systems known to exhibit non-KAM 
chaos [13, EH- Its rich dynamical properties have been 
studied by many authors and provide insights for un- 
derstanding a wide range of problems in, for example, 
plasma physics, tokamak fusion, turbulent fluid dynam- 
ics, ion traps and quasi-crystals [3, [l?], EH- Some in- 
teresting effects of dissipation on non-KAM chaos, in- 
cluding the identification of Cantor sets in phase space, 
were studied in Ref. [l||. But, to our knowledge, the 
bifurcation phenomena considered in the present paper, 
and their relation to resonance effects in SLs, have not 
been considered elsewhere. Moreover, what has not been 
clear previously is the nature of the characteristic insta- 
bilities and dynamical regimes that can be induced by 
dissipation in systems that exhibit non-KAM chaos, or 
how these regimes evolve with variation of the control 
parameters. Our analysis of such dynamics for electrons 
in a SL with a tilted magnetic field is therefore also of 
general interest in nonlinear dynamics and in the diverse 
areas of physics mentioned above, which involve the be- 
havior of a harmonic oscillator driven by a plane wave. 



III. EFFECT OF DISSIPATION ON CHAOTIC 
DYNAMICS. 

In Ref. [l3T] . the dynamics of a miniband electron in an 
electric and tilted magnetic field were studied in the ab- 
sence of dissipation by solving Eq. ^ with a = 0. Two 
distinct types of chaotic trajectories were identified from 



the different patterns that they produce in the Poincare 
sections through phase space. The first type of pattern 
occurs only when the resonance condition 

uj b = ruj c (9) 

is satisfied, where r is a rational number. At such reso- 
nances, the chaotic orbits map out a stochastic web com- 
prising an intricate mesh of filaments, which threads the 
phase space and is a unique feature of non-KAM chaos. 
By contrast, the second type of pattern, a large continu- 
ous "chaotic sea" in the Poincare section, does not require 
the resonance condition to be satisfied. 

We now illustrate these two types of chaos and inves- 
tigate how the electron dynamics and associated phase 
space patterns change when dissipation is introduced by 
increasing a from zero in Eq. (JSJ). First, we consider the 
case when the phase space of the dissipation-free (a = 0) 
system contains a stochastic web corresponding to the 
r = 1 resonance condition attained when F = 2.9 x 10 a 
Vm" 1 , B = 2 T, and 6 = 30°. The black dots in Fig. 
[2fa), which merge to form a continuous pattern, show a 
stroboscopic Poincare section constructed from the elec- 
tron trajectories by plotting the momentum components 
iPyiPz) in the plane of the SL layers at discrete times sep- 
arated by the Bloch period Tb — 2tt/lub- The Poincare 
section reveals a stochastic web, formed by the black dots 
and labeled "0" in Fig. Ufa), which contains both ring- 
like and radial filaments. The radial filaments act as 
conduction channels, which enable the electron to diffuse 
rapidly outwards away from the web centre, thus increas- 
ing its momentum pl = (Py + P 2 ) 1 ^ 2 in the plane of the 
SL layers. In real space, the electron progresses rapidly 
along the SL axis [241 ] . gaining kinetic energy from the 
electric field, which is transferred into the y — z motion 
by the tilted magnetic field [HI, [3]. To illustrate the 
resonant derealization of the electron trajectories that 
results from stochastic web formation, the black curve 
labeled "0" in Fig. [2Kb) shows an electron orbit within 
the stochastic web projected onto the x~z plane in which 
the magnetic field lies (axes inset). The electron starts 
from rest at the left-hand edge of the orbit, which corre- 
sponds the center point (p y ,p z ) — (0,0) of the stochastic 
web in Fig. [2ja). A slow modulation of the density of or- 
bital loops, which gradually increases with increasing x, 
is the only indication of irregularity in the electron tra- 
jectory labeled 0. For this reason, stochastic web chaos is 
often called "weak chaos" [l8[ , because the narrow width 
of the web filaments largely suppresses the highly erratic 
behavior usually found for nonintegrable systems. As 
the electron travels along the orbit from left to right in 
real space, it diffuses outwards along the almost vertical 
radial filament in the stochastic web towards higher p z 
values. The web filaments enmesh islands of stability, 
which are slices through invariant phase space tori gen- 
erated by regular quasiperiodic trajectories [HI, [HI- For 
clarity, in Fig. [2^a) we show no tori corresponding to 
orbits within these islands of stability. 

The inclusion of even a very small dissipation term in 
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FIG. 2: (Color) (a) Black dots: stroboscopic Poincare section 
calculated for the chaotic trajectory executed by an electron 
starting from rest in the dissipationless limit a — 0. This sec- 
tion reveals the inner filaments of a stochastic web. Colored 
triangles numbered 1-6: Poincare sections calculated for six 
coexisting stable limit cycles formed when a — 10 10 s^ 1 . Blue 
circle numbered 7: Poincare section calculated for the single 
stable limit cycle found when a — 10 12 s _1 . (b) Electron tra- 
jectories corresponding, by number and color, to the labeled 
features in the Poincare sections shown in (a). All orbits are 
projected onto the x — z plane (axes inset), have a common 
spatial scale shown by the horizontal line whose length is 30 
SL periods, and are plotted over a fixed time interval of 9.45 
ps. The field parameters B = 2 T, 9 = 30° and F = 2.9 x 10 5 
Vm" 1 satisfy the resonance condition lob/u)c = L 



Eq. © destroys both the stochastic webs and islands 
of stability, replacing them with attracting limit cycles 
to which all orbits tend as time progresses. To illustrate 
this, we first set a = 10 10 s _1 , which is far less than the 
momentum relaxation rate W m ~ 10 12 — 10 13 s _1 typical 
of electrons in semiconductor SL devices [l3l [T3 |. For 



this value of a, there are several distinct limit cycles, six 
of which are marked by the numbered and colored trian- 
gles in the stroboscopic Poincare section shown in Fig. 
UKa). Each limit cycle appears within, and attracts tra- 
jectories from, one of the islands of stability found when 
a = 0. Segments of the electron trajectories for each 
limit cycle are shown in Fig. [2Kb) . Each orbit is shown 
over the same fixed time interval of 9.45 ps, chosen when 
t is large enough for the orbit to approach very close to 
the limit cycle, and linked, by number and color, to the 
corresponding points in the phase space plot shown in 
Fig. [Ha). The orbital segments reveal that the inclu- 
sion of weak dissipation makes the electron trajectories 
more regular because each limit cycle creates a particular 
repeating loop pattern. 

As the level of dissipation increases, the different at- 
tracting limit cycles degenerate until only one remains 
in the phase space, whose position for a — 10 12 s _1 is 
marked by the blue circle labeled "7" in Fig. [Ha). Part 
of the corresponding electron orbit is shown by the blue 
curve labeled "7" at the bottom of Fig. (U^b). The sim- 
ple loop structure of this orbit transports the electron 
rapidly through the SL and is qualitatively similar to 
that of orbit 0, which generates the stochastic web in 
the absence of dissipation. So, although dissipation com- 
pletely changes the phase space structure, replacing the 
stochastic web with a small number of limit cycles, it 
preserves the overall form of the delocalized orbits found 
when ljb = rCo c and actually enhances resonant electron 
transport by preventing the spatial compression of the 
orbital loops, which slows the electron in the absence of 
dissipation (see, for example, orbit in Fig. [2Kb)). 

We now consider the second type of Hamiltonian chaos, 
which, in the absence of dissipation, produces an ex- 
tended chaotic sea in phase space even when lob/OJc is 
not a rational number. The Poincare section formed by 
the black dots in Fig. [3Ka) shows the chaotic sea gener- 
ated by solving Eq. © for an electron starting from rest 
in the SL with a = 0, F = 7.5 x 10 5 Vm" 1 , B = 4.75 T, 
and 9 — 45°, for which r « 1.33. Introducing dissipation 
into the equation of motion affects this type of extended 
chaotic sea far less than stochastic webs. In particular, 
although dissipation changes the form of the chaotic sea 
and the orbits within it, the orbits remain unstable and 
spatially irregular even for fairly large values of a. The 
reason is that dissipation transforms the chaotic sea into 
a chaotic attractor, rather than into a stable limit cycle 
like those labeled 1-6 in Fig. 2. In Fig. O^b), the yellow, 
red and green dots show the chaotic attractors generated 
by plotting (p y ,p z ) at every Bloch period for electron tra- 
jectories starting from rest with, respectively, a = 10 10 
s _1 , 10 11 s _1 , and 10 12 s _1 . The irregular arrangement 
of points within these attractors demonstrates that they 
represent chaotic orbits, a conclusion that is confirmed 
by calculating their Lyapunov exponents. However, with 
increasing a, the dimension of the chaotic attractors di- 
minishes. For the three attractors shown in Fig[3Kb), the 
Lyapunov dimension calculated using the Kaplan- Yorkc 
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FIG. 3: (Color) (a) Black dots: stroboscopic Poincare sec- 
tion calculated for an electron trajectory starting from rest 
in the dissipationless limit a — 0. Colored dots: stroboscopic 
Poincare sections calculated for the chaotic attractors of elec- 
trons starting from rest when a = 10 10 s -1 (yellow dots), 
10 11 s" 1 , (red dots), and 10 12 s" 1 (green dots), (b) Electron 
trajectories corresponding, by specified value of a and color, 
to the Poincare sections shown in (a) . All orbits are projected 
onto the x — z plane (axes inset), have a common spatial scale 
shown by the horizontal line whose length is 10 SL periods, 
and are plotted over a fixed time interval of 22.68 ps, which 
is long enough to reveal the form of the electron trajectories. 
B = 4.75 T, 9 = 45°, F = 7.5 x 10 5 Vm" 1 , r « 1.33. 

formula [H,[26| is 2.987, 2.887, and 2.386. Chaotic elec- 
tron orbits calculated for a = (black curve), 10 10 s _1 
(yellow curve), 10 11 s _1 (red curve), and 10 12 s _1 (green 
curve) are shown projected into the x — z plane in Fig. 
(3fb). Comparison of the trajectories and Poincare sec- 
tions shown in Fig. [3] (a) and (b) for different values of 
a reveals that as the dissipation increases, the chaotic 
attractors gradually contract towards simple curves in 
phase space and, as a consequence, the orbits become 
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spatially more regular. We emphasize that, in spite of 
the presence of energy dissipation, the electron trajec- 
tories remain chaotic, although their their topological 
structure in phase space changes as a varies over a wide 
range of values. In contrast to the Hamiltonian case, the 
phase space of the dissipative system contains unstable 
limit sets, which attract trajectories from certain regions 
of phase space known as basins of attraction. However, 
for very large a > 10 13 s _1 chaos completely disappears 
leaving only periodic orbits in the system. 



IV. STABILITY AND BIFURCATION 
ANALYSIS. 

In this section, we investigate how the stability of the 
electron orbits changes with the field parameters F, B, 
and 9, and explore the bifurcation sequences that drive 
the transition between regular motion and dissipative de- 
terministic chaos. We consider a fixed value of a = 10 12 
s _1 . First, we identify different regimes of electron dy- 
namics by calculating the Lyapunov exponents for those 
limit sets that attract trajectories starting from rest. To 
characterize the stability of the electron motion, we cal- 
culate the largest Lyapunov exponent, A. The orbit is 
chaotic if A > 0, but otherwise is stable and regular. 

The color map in Fig. 0Ja) shows the value of A cal- 
culated as a function of F and B for — 45°. For small 
F or B, the system exhibits stable periodic oscillations, 
of frequency u>b, driven by the plane wave. We call such 
oscillations "Period-1 limit cycles". As F or B grows, the 
system undergoes a sequence of bifurcations, which even- 
tually leads to the appearance of chaos within the yellow 
and red areas in Fig. 0Ja), where A > 0. The crosses in 
Fig. HJa-c) mark the coordinate in F — B space for which 
an electron starting from rest eventually approaches the 
chaotic attractor shown green in Fig(3Jb). The location 
and shape of the various islands of stable and chaotic elec- 
tron motion in the F — B plane [Fig. [UJa)] can be under- 
stood by considering the bifurcations that the limit cycles 
undergo as F and/or B increases. The black curves in 
Fig. 0Jb) show the loci of the primary sequence of period- 
doubling bifurcations, i.e. those that occur first as F in- 
creases. Each locus encloses a large "V" -shaped region, 
which we call a "period-doubling tongue" surrounding 
at least one area of chaos in the F — B plane. When 
the field parameters enter one of these period-doubling 
tongues, the limit cycle electron trajectory undergoes a 
period-doubling bifurcation. As a result of this bifurca- 
tion, the initial limit cycle loses its stability and, in its 
vicinity, another stable limit cycle appears whose period 
is twice as long. 

A second type of bifurcation, namely the saddle-node 
bifurcation for limit cycles, also plays a key role in the 
electron dynamics. Loci of the primary saddle-node bi- 
furcations are shown by the green curves in Fig. [U[b), 
which enclose "saddle-node" tongues indicated by the 
green shaded areas. The top edges of these tongues 
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FIG. 4: (Color) (a) Color map showing A calculated as a 
function of F and B for the attractors approached by elec- 
tron trajectories starting from rest. The color scale is shown 
on the right, (b) Black (green) curves are loci of the primary 
period-doubling (saddle-node) bifurcations of limit cycles cor- 
responding to resonances with rational r values (marked for 
r — 1,1/2,1/3). Light green shaded areas are saddle-node 
tongues, (c) Location of chaotic attractors (red areas where 
A > 0) and additional bifurcations within the region of F — B 
space bounded by the r — 1 primary period-doubling bifur- 
cation (outer black curve). Black (green) curves are loci of 
period-doubling and saddle-node bifurcations corresponding 
to rational values of r ^ 1. Light green shaded areas are 
saddle-node tongues. Dashed line marks B — 2 T, and blue 
arrow shows the path through F — B space corresponding to 
the bifurcation diagram in Fig. 5. In all panels, crosses mark 
(F, B) coordinate corresponding to the chaotic orbit shown 
(green) in Fig. 3(b). 9 = 45° and a = 10 12 s" 1 . 



are close to those of the period-doubling bifurcation 
loci [black curves in Fig. BJb)] and lie near the lines 
F = r(eHcos6/m*ed)B (r = 1, 1/2, 1/3, ...) along which 
lub = tu) c . Consequently, each of the period-doubling 
and saddle-node tongues can be associated with a par- 
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ticular value of r [those for r = 1, 1/2, 1/3, ... are labeled 
in Fig. E(b)], which relates to a resonance of a Period-1 
cycle. Note, however, that when a ^ 0, nonlinear res- 
onances can occur near, but not exactly at, field values 
for which r is a rational number. When the field coor- 
dinate (F, B) enters a saddle-node tongue, a new stable 
limit cycle is created, which corresponds to an unbounded 
electron trajectory, like that labeled "7" in Fig. [2(a) for 
r = 1. This new limit cycle is twinned with an unstable 
saddle counterpart. As the field coordinate (F, B) leaves 
the boundary of the saddle-node tongue, the stable and 
saddle cycles approach one another in phase space and 
eventually combine via the saddle-node bifurcation that 
occurs at the edge of the tongue (we consider this process 
in more detail below). At this bifurcation point, the sta- 
ble and unstable saddle cycles both disappear, and the 
electron trajectory jumps abruptly onto another stable 
orbit located in a different part of phase space. 

Within the area of F — B space bounded by each of 
the primary period-doubling tongues in Fig. 3(b), the 
system undergoes a cascade of additional bifurcations. 
In Fig. [4(c), we illustrate this for the region of F — B 
space enclosed by the r — 1 primary period-doubling 
bifurcation (outer black curve). Within this region, fur- 
ther period-doubling and saddle-node bifurcations occur 
along the black and green curves respectively. For exam- 
ple, the sharp saddle-node and period-doubling tongues 
whose apexes are at (F, B) « (5 x 10 5 Vm _1 , 2 T) in Fig. 
[4(c) are associated with the r = 2 nonlinear resonance 
and produce orbits of period 2Tb- Due to the presence 
of multiple tongues, the bifurcation diagram exhibits a 
complex self-similar structure. It is characterized by a 
cascade of bifurcations, which drives into chaos the elec- 
tron trajectories that start from rest. These chaotic tra- 
jectories lie within the yellow and red areas of Fig. 0(a) 
and (c), where A > 0. 

An interesting effect occurs within those regions of 
F — B space where period-doubling and saddle-node 
tongues intersect, for example the area marked by the 
cross in Fig. 0(c). As noted above, each period-doubling 
tongue contains a family of attractors. But in addition, 
the corresponding saddle-node tongue contains a distinct 
stable limit cycle and a saddle (unstable) periodic orbit. 
Consequently, in the overlap region, at least two differ- 
ent attracting limit sets can coexist, which means that 
multistability occurs in the system. The initial condi- 
tion of the orbit determines which attractor it eventually 
approaches. 

To reveal in detail the bifurcations that the system un- 
dergoes inside the tongues in Fig. 0Jc), we now consider 
how the electron orbits change as B increases along the 
blue arrow in Fig. 2(c) with fixed F = 3x 10 5 V m -1 . In 
particular, we construct the single-parameter bifurcation 
diagram shown in Fig. El by finding the limit sets that at- 
tract electrons starting from rest for 0.5 < B < 3 T and, 
for each B, plotting the values that p z attains (black 
dots in the figure) whenever p z = in the limit cycle. 
For B < 1.15 T, this attractor is a Period-1 cycle, which 
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FIG. 5: (Color) (a) Poincare sections constructed by plotting 
p z whenever p z — along limit cycles for electrons starting 
from rest (black dots) and for stable/unstable Period-1 limit 
cycles (red/green circles) with F = 3x 10 5 Vm _1 , = 45° and 
a — 10 12 s _1 . Increasing B induces period-doubling (PD1, 
PD2) and saddle-node (SN1, SN2) bifurcations, which corre- 
spond to crossing the black and green curves along the blue 
arrow in Fig. 4(c), and a boundary crisis (BC). Arrow marks 
B value for which r = 1. (b) An enlargement of the Poincare 
section in the vicinity of saddle-node bifurcation SN1 and 
boundary crisis BC. 



produces a single black dot in Fig. [5(a) at a p z value that 
gradually increases with increasing B. When B ~ 1.15 
T, the Period-1 cycle undergoes a period-doubling bifur- 
cation, labeled PD1 in Fig. [5(a), which corresponds to 
entering the r = 1 primary period-doubling tongue en- 
closed by the outer black curve in Fig. BJc). As a result, 
p z takes one of two values when p z — 0, causing the black 
dots to split into two distinct branches in Fig. [5(a). As 
B increases further, the limit cycle for electrons start- 
ing from rest undergoes a sequence of additional period- 
doubling bifurcations, leading to the onset of chaos when 
B ~ 1.85 T. Thereafter, chaotic motion alternates with 
"stability windows" , which is a characteristic feature of 
chaos induced by a cascade of period-doubling bifurca- 
tions [26[, and suddenly disappears when B rj 2.461 T, 
near the r=l nonlinear resonance, which occurs at the B 
value marked by the arrow in Fig. E(a). This suppres- 
sion of chaos occurs because the attractor for electrons 
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starting from rest switches back to a stable Period- 1 cy- 
cle owing to a "boundary crisis" , which we explain in 
detail below. At this transition from chaotic to stable 
motion, the distribution of black dots in Fig. E{a), which 
is shown more clearly in the enlargement in Fig. [5fb), 
changes from an irregular scatter to an ordered set of 
points located along a single-valued curve. 

To gain further insights into the electron dynamics, we 
now consider the nature and evolution of Period- 1 cycles 
in the system over the entire range of B values shown in 
Fig. OJa). To do this, we plot colored circles in Fig. [5] 
showing the values of p z whenever p z = in the Period- 1 
cycle. If the Period- 1 cycle is stable (unstable) we plot a 
red (green) circle. For B < 1.15 T, the Period-1 orbit is 
stable and is the only attractor in the system, being the 
limit cycle that all trajectories approach irrespective of 
their initial conditions. 

At the first period-doubling bifurcation [labeled PD1 
in Fig. [S{a)], the Period-1 cycle becomes unstable (and 
therefore no longer an attractor) and co-exists with the 
period-doubled stable limit cycle produced by the bi- 
furcation. When B reaches w 2.45 T [labelled SN1 in 
Fig. E{a) and (b)] it enters the r — 1 primary saddle- 
node tongue shown in Fig. HJc). The associated saddle- 
node bifurcation creates new stable and unstable saddle 
Period-1 cycles, which generate, respectively, the red and 
upper green circles in Fig. EJa),(b). 

Thus, the two new Period-1 cycles formed by the 
saddle-node bifurcation coexist with the original Period- 
1 cycle, which produces the lower colored circles in Fig. 
[5j This original Period-1 cycle stabilizes (circles change 
from green to red) when B reaches ~ 2.67 T and there- 
fore crosses the upper edge of the r = 1 primary period- 
doubling tongue shown in Fig. IDJc). At a slightly higher 
B value of « 2.73 T, it coalesces with the unstable sad- 
dle Period-1 cycle and thereby disappears as a result of 
saddle-node bifurcation SN2. At this bifurcation, the 
field coordinate (F, B) exits through the upper edge of 
the r = 1 primary saddle-node tongue in Fig. IHc). 

At higher B values there is a single Period-1 cycle, 
which is the only attractor in the system. Note that 
the nonlinear r = 1 resonance has a pronounced effect 
on the bifurcation diagram in Fig. [5] In particular, it 
gives rise to the characteristic loop formed by the dif- 
ferent branches of the Period-1 limit cycles that coexist 
between bifurcations SN1 and SN2 in Fig. [5ja) and also 
causes the abrupt disappearance of the chaotic attractor 
at the boundary crisis, BC. We now consider the latter 
in more detail. 

As explained above, Fig. [5] reveals a small field range 
between B « 2.450 T and B w 2.461 T (i.e. between 
the points labeled SN1 and BC in the figure) where 
two distinct attractors coexist, a chaotic one and the 
stable Period-1 limit cycle. Figure [l^a) shows a stro- 
boscope Poincare section calculated for the chaotic at- 
tractor (black dots), stable Period-1 limit cycle (red cir- 
cle), and unstable saddle Period-1 limit cycle (green cir- 
cle) constructed by plotting the momentum components 
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FIG. 6: (Color) Stroboscopic Poincare sections calculated at 
times t = IT B (I = 0, 1,2, ...) when B « (a) 2.46 T just be- 
fore boundary crisis, (b) 2.461 T at the boundary crisis. The 
sections are constructed by plotting (p y ,p z ) whenever p z — 
along the chaotic attractor for electrons starting from rest 
(black dots), and for the stable Period-1 limit cycle (red cir- 
cles) and unstable saddle Period-1 limit cycle (green circle), 
which are formed by saddle-node bifurcation SN1 in Fig. 5. 
White (gray) regions are basins of attraction for the chaotic 
attractor and stable Period-1 limit cycle respectively. Ar- 
rows in (a) pointing towards (away from) the green circle 
represent stable (unstable) manifolds of the unstable saddle 
Period-1 limit cycle. Arrow and dotted curve in (b) show 
the direction along which the electron trajectory eventually 
leaves the chaotic attractor and subsequently approaches the 
stable Period-1 limit cycle. F = 3 x 10 5 Vm _1 , 6 = 45° and 
a = 10 12 s- 1 . 



(p v ,p z ) at times t = IT B (I = 0, 1,2, ...) when B = 2.46 
T, i.e. just before the boundary crisis that destroys the 
chaotic attractor. Each attractor has a distinct basin 
of attraction comprising the set of initial conditions for 
which the trajectory will eventually approach the attrac- 
tor. The basins of attraction for the chaotic set and 
Period-1 limit cycles are shown respectively by the white 
and gray areas of the p y — p z plane in Fig. [6ja). Note 
that the boundary between these two basins of attraction 
is the stable manifold of the unstable saddle periodic or- 
bit (green circle), which originates from the saddle-node 
bifurcation labeled SN1 in Fig. [5] Trajectories with ini- 
tial conditions exactly on this stable manifold approach 
the unstable saddle periodic orbit along the directions 
shown by the arrows pointing towards the green circle in 
Fig. [6t[a) . Conversely, the unstable manifold of the sad- 
dle, represented schematically by the two arrows pointing 
away from the green circle in Fig. [|3a) repels electron tra- 
jectories away from the vicinity of the saddle orbit. Note 
that the lower right-hand part of the chaotic attractor in 
Fig. EUa) approaches very close to the boundary of its 
(white) basin of attraction. This suggests that if chang- 
ing B increases the size of the chaotic attractor, there will 
be a boundary crisis, which occurs when an attractor hits 
the boundary of its basin of attraction. 

To investigate this possibility, we studied how the 
phase space structure shown in Fig. [Sfa) responds to 
small changes of B. Our calculations confirm that in- 



9 



creasing B to w 2.461 T does indeed trigger a bound- 
ary crisis, which causes the abrupt disappearance of the 
chaotic attractor at the field value labeled BC in Fig. 
[5j To illustrate this, Fig. Hfb) shows the stroboscopic 
Poincare section calculated for the chaotic set (black 
dots), stable Period- 1 limit cycle (red circle) and unstable 
saddle Period- 1 limit cycle (green circle) at the boundary 
crisis. The lower right-hand edge of the chaotic attractor 
touches the saddle (green circle), which, for lower B, was 
located on the boundary between the basins of attrac- 
tion of the chaotic and Period-1 attractors [Fig. [DJa)]. 
When the chaotic attractor and saddle point touch at 
the boundary crisis, all trajectories that approach the 
chaotic attractor eventually leave it at the bottom right- 
hand tip in Fig. [6]Jb) , then cross the curve in the p y — p z 
plane that used to separate the two basins of attraction, 
and thereafter move in the direction of the arrow along 
the dotted curve in Fig. [H^b) towards the stable Pcriod-1 
limit cycle (red circle). Consequently, for all initial con- 
ditions, the electron trajectory eventually approaches the 
stable Period-1 limit cycle, either directly or via the tip of 
the former chaotic attractor, which now allows the elec- 
tron trajectories to escape. Over the field range between 
the boundary crisis and the second period-doubling tran- 
sition (i.e. between the fields labeled BC and PD2 in Fig. 
[5]) the Period-1 limit cycle is therefore the only attractor 
in the system and its basin of attraction [gray in Fig. 
E]Jb)] fills the entire phase space. 

We now consider how the attractors evolve as B sweeps 
either up or down through the boundary crisis. As noted 
above, for B values in the range « 1.15 to sa 2.45 T (i.e. 
between points PD1 and SN1 in Fig. [5]), the chaotic set 
(black dots in Fig. [5]) is the only attractor in the system. 
As B increases above w 2.467 T, the state of the sys- 
tem remains pinned to the chaotic attractor even though 
a second attractor (red circles in upper part of Fig. [5]) 
emerges from the saddle-node bifurcation labeled SN1 in 
Fig. [5] When B reaches the value at which the bound- 
ary crisis occurs, the system state changes abruptly from 
chaotic to the Period-1 attractor, as discussed above. 
Conversely, if B now decreases, the state remains pinned 
to the Period-1 attractor until the saddle-node bifurca- 
tion SN1 [where the red and green circles meet in Fig. 
[SJb)] removes this attractor when B m 2.45 T, i.e. at 
a field value below that corresponding to the boundary 
crisis. Consequently, the transition between chaotic at- 
tractors and regular limit cycles occurs at a higher B 
value when the field sweeps up through the boundary 
crisis than when it sweeps down, and the system exhibits 
hysteresis. For given F, the hysteresis occurs over the 
magnetic field range for which the r = 1 saddle-node 
tongue, shown by the upper gray region in Fig. HJc), 
overlaps with the upper red island in the figure, in which 
the attractor is chaotic. Figure|3Jc) shows that this over- 
lap region broadens with increasing F, meaning that the 
hysteresis becomes more pronounced, that is, it extends 
over a wider magnetic field range. The occurrence of hys- 
teresis in bifurcation diagrams like Fig. [5] is a common 



feature of resonances in non-linear systems [27] . 

Note that two distinct mechanisms cause the onset of 
chaotic limit cycles in the system. As B increases from 0, 
the transition to chaos occurs via the cascade of period- 
doubling bifurcations shown in Fig. [5J Conversely, as B 
decreases from 3 T, chaos switches on abruptly at the 
saddle-node bifurcation (point SN1 in Fig. [5} that re- 
moves the stable Period-1 limit cycle. 
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FIG. 7: (Color) (a) Color map showing A calculated as a func- 
tion of F and 9 for the attractors approached by electron tra- 
jectories starting from rest when B = 2 T and a — 10 12 s _1 . 
Color scale is shown on the right. Yellow (green) curves are 
loci of period-doubling (saddle-node) bifurcations of electron 
limit cycles. Cross marks (F, 6) coordinate corresponding to 
limit cycle 7 in Fig. 2. 



The two main mechanisms responsible for the onset of 
chaos (the period-doubling sequence and saddle-node bi- 
furcation) occur for a wide range of 9 values. This can 
be seen from Fig. \7\ which shows the loci of the main 
period-doubling (yellow curves) and saddle-node (green 
curves) bifurcations in the F — 9 plane calculated for fixed 
B = 2T. These loci are overlaid on a color map show- 
ing the value of the largest Lyapunov exponent for the 
limit cycles approached by electron trajectories starting 
from rest. Note that, as in Fig. [U the limit cycles be- 
come chaotic (yellow or red in the color map) within the 
regions of F — 9 space bounded by the period-doubling 
bifurcations (yellow lines), which partially overlap the ar- 
eas enclosed by the saddle-node bifurcations (green lines). 
As a consequence of this overlap, multistability, and the 
associated hysteretic evolution of limit cycles, occurs for 
a range of F and 9 in the same way that hysteresis oc- 
curs within islands in the F — B plane (see Fig. 2]). Note, 
however, that for the particular field value B = 2 T used 
to calculate Fig. chaos only occurs for 9 > 40°. For 
smaller values of 9, electrons starting from rest approach 
only stable limit cycles, for example that labeled 7 in Fig. 
2, which is calculated for the (F, 9) coordinate marked by 
the cross in Fig. [Tj 
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FIG. 8: (Color) (a) Amplitude, A, of p z (t) oscillations cal- 
culated versus F for a stable/unstable Period- 1 limit cycle 
(red/green circles) and a stable/unstable Period-2 limit cycle 
(red/green triangles). Arrows mark positions of r = 1 and 
2 resonances, (b) (v x ) calculated as a function of F for the 
attractor approached by electrons starting from rest. Peaks 
labeled R1(R2) occur at F values near the r = 1(2) nonlinear 
resonances. Labels 01(02) mark off-resonance field values 
for which r = 1.52(2.53). (c) Orbits approached by electrons 
starting from rest at the resonant (R1,R2) and off-resonant 
(01,02) F values marked in (b). All attractors are projected 
onto the x — z plane (axes inset) , have a common spatial scale 
shown by the horizontal line whose length is 30 SL periods, 
and are plotted over the same time interval. B — 2T, 9 = 45° , 
a = 10 12 b- 1 . 



V. THE EFFECT OF RESONANCES ON 
ELECTRON VELOCITY 

In Section II, we showed that electron motion in the SL 
can be described by Eq. (J5|) , corresponding to a damped 



harmonic oscillator driven by a plane wave. When the 
driving frequency, lob oc F, is commensurate with the 
natural frequency, uJ c oc B, of the harmonic oscillator 
nonlinear resonances occur in the system. The effect of 
such resonances on the electron dynamics is complex. For 
example, in Sections III and IV we saw that the r = 1 res- 
onance creates a stable limit cycle, which corresponds to 
a spatially unbounded electron trajectory, via a saddle- 
node bifurcation, and also triggers a boundary crisis. In 
this section, we investigate how resonances affect the mo- 
mentum and drift velocity of the electrons, both parallel 
and perpendicular to the SL layers. We focus on these 
dynamical variables because they strongly influence the 
transport characteristics, such as I(V) curves, measured 
in experiments (l4| . 

First, we consider how resonances affect the variation 
of the in-plane momentum component, p z , as the elec- 
tron performs a limit cycle. This momentum compo- 
nent is the most natural parameter to study because it 
is the dependent variable in Eq. |6j). For each limit 
cycle, we define the amplitude of the p z oscillations to 
be A = (pf ax -p™ m )/2, where pf ln and pf ax are, re- 
spectively, the minimum and maximum values of p z at- 
tained during the limit cycle. We then investigate how A 
varies with F for several different limit cycles, including 
those produced by the bifurcations induced by changing 
F, In Fig. [5Ja), we show A calculated as a function of F 
for Period- 1 (circles) and Period-2 (triangles) limit cycles 
when B = 2 T and 6 = 45°. The A(F) plots, analogous 
to the frequency response curves familiar to engineers, 
exhibit the resonant features that characterize nonlinear 
dissipative systems. The resonances are clearl y s een as 
sharp increases of A when r is close to rational [28[ , with 
the largest peaks [arrowed in Fig. Ufa)] occurring near 
the F values for which r = 1 and 2. The smaller peaks 
seen at lower F correspond to rational r < 1. A Period- 1 
limit cycle exists over the entire range of F shown in 
Fig. EJa), but undergoes bifurcations as F increases, 
which changes it between stable (red circles) and unstable 
(green circles) . The location of the primary bifurcations 
for B = 2 T and 6 = 45° can be seen from Fig. [3 For F 
values close to 2.35 x 10 5 Vm -1 , when r ~ 1, there is a 
saddle-node tongue (enclosed by the green curve in Fig. 
[7]), which means that up to three Period- 1 orbits coexist 
(as in the loop region of Fig. [5]). Consequently, the A(F) 
curve for the Period- 1 limit cycle contains a loop when 
r«l. Within the loop region, the stable Period- 1 limit 
cycle undergoes a period-doubling bifurcation (along the 
right-hand yellow curve in Fig. [7J, which produces the 
Period-2 limit cycle whose A(F) curve [triangles in Fig. 
IHa)] reveals large resonant peaks when r = 1 and 2. 

The electron velocity component, v x , along the SL axis 
is a particularly important dynamical variable because 
it determines the electrical current that flows through 
the device and the power dissipation rate per particle, 
Pd = eFv x . It follows from the equations of motion (J5]) 
and ([5]) that for each attractor, v x oscillates as a function 
of time. To determine the drift component of the elec- 
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FIG. 9: Normalized Ii(v x ) curves (with gray underfill) calcu- 
lated for electron trajectories labeled in Fig. [8] (b) and (c) as: 
(a) Rl, (b) Ol, (c) R2, (d) 02. Vertical dashed lines mark 
extremal velocities v x = ±%«i. B — 2 T, — 45°, a = 10 12 
s" 1 . 



tron velocity, which determines the current flow [1 41 ] . we 
therefore calculate the time average (v x ) of the velocity 
throughout the attractor. Figure HJb) shows (v x ) calcu- 
lated as a function of F for the attractor approached by 
electrons starting from rest with B — 2 T and 6 = 45°. 
The shape of this graph is similar to the A(F) data shown 
in Fig. EI a). In particular, there is strong resonant en- 
hancement of (v x ) near rational r values where maxima 
also occur in the A(F) curves calculated for the Period-1 
and Period-2 limit cycles. The peak labeled Rl in the 
(v x ) versus F curve corresponds to the r = 1 resonance 
of the Period-1 limit cycle, whereas peak R2 originates 
from the r = 2 resonance of the Period-2 limit cycle. 
Those peaks to the left of Rl arise from resonances of 
the Period-1 limit cycle. To relate the resonant peaks in 
the (v x ) versus F curve directly to the electron motion, 
Fig. Etc) shows segments of Period-1 (lower and two up- 
per trajectories) and Period-2 (third orbit down) limit 
cycles calculated at the field values labeled Rl, R2, Ol, 
and 02 in Fig. E^b) over a fixed time interval of 12.38 
ps. During this time interval, the electron travels much 
further along in the two orbits on resonance (Rl and R2) 
than in the two off-resonance orbits (Ol and 02). 

To further investigate how resonances affect the elec- 
tron transport, we calculated the probability distribu- 
tion, H(v x ), of all the v x values attained during one pe- 
riod of the limit cycles approached by electrons starting 
from rest both on and off resonance. Figure [9] shows 
normalized H(v x ) curves calculated for the orbits labeled 
(a) Rl, (b) Ol, (c) R2, (d) 02 in Fig. EI (b) and (c). All 
attainable v x values lie between the two vertical dashed 
lines located at v x = ±v max , where v max = dA/2h is the 
maximum velocity of an electron in the lowest miniband. 
Comparison of the velocity distributions on (a),(c) and 
off (b),(d) resonance reveals that on resonance, the most 
probable v x value approaches v max , thus producing the 



FIG. 10: (v x ) calculated as a function of increasing (black 
curve) and decreasing (under-filled gray curve) B (upper 
axis). Lower axis shows values of r a 1/B. Peaks in (v x ) 
originate from r = 1,3/2, and 2 resonances (located by ar- 
rows). The difference between the two curves originates from 
the hysteretic evolution of the electron limit cycle discussed 
in Section IV. Dashed curve shows Vd (with no vertical offset) 
calculated as a function of B (or r) using Eq. (|10[) . F = 10 6 
Vm-\ 6 = 45°, a = 10 12 s" 1 . 



large peak (v x ) values shown in Fig. [5Jb). 

The hysteretic response exhibited by the limit cycles 
as B is swept up and down through a boundary crisis 
has a clear manifestation in the (v x ) versus B curves 
calculated for large F values where the hysteresis is most 
pronounced, as discussed in Section IV. To illustrate this, 
Fig. [TUl shows (v x ) calculated as a function of increasing 
(black curve) and decreasing (under-filled gray curve) B 
(upper axis) when F = 10 6 Vm" 1 and = 45°. Note that 
the lower axis of the graph shows the ratio, r oc 1/B, of 
the Bloch and cyclotron frequencies. Near the r = 1, 3/2, 
and 2 resonances (marked by arrows in Fig. ITTTj) . saddle- 
node tongues coexist with chaotic attractors. This can 
be seen, for example, from Fig. 0|c), which, as discussed 
in Section IV, reveals a large overlap between the r = 1 
saddle-node tongue and upper red island of chaos when 
F = 10 6 Vm -1 . As a consequence, the limit cycle evolves 
in a different way when B increases or decreases, and so 
the corresponding (v x ) versus B curves exhibit hysteresis. 
This hysteresis is particularly pronounced near the r = 1 
and r = 2 resonances in Fig. [101 which reveals that for 
some B values (v x ) differs by a factor of two between up 
and down sweeps of B. Since the current flow depends 
critically on the electron drift velocity, we predict that 
measured 1(B) curves should also reveal strong resonant 
peaks that exhibit clear hysteresis. 

In Fig. 1101 we compare the (v x ) versus B curves ob- 
tained by including dissipation a priori in the equations 
of motion, with drift velocities, Vd{B\ calculated using 
the traditional Esaki-Tsu approach [3j, based on non- 
dissipative electron trajectories. Within this approach, 
which we used in our previous studies of electron trans- 
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where = x is determined from Eqs. ([6]) and (8), set- 
ting a = 0, and the electron scattering time, r, includes 
contributions from both elastic and inelastic scattering 
processes [14| . The Vd(B) variation obtained from Eq. 
(TTU|) is shown by the dashed curve in Fig. [TOJ Com- 
parison with the {v x } versus B curves (black and under- 
filled gray) reveals that both approaches give qualita- 
tively similar results, with prominent maxima at B val- 
ues close to r = 1, 3/2, and 2. However, the positions 
of these maxima do not always coincide exactly, in par- 
ticular for r = 1, mainly because including dissipation a 
priori in the equations of motion shifts the resonances 
slightly away from rational values of r, as discussed in 
Section IV. Although the resonant peaks in the (v x ) ver- 
sus B and Vd(B) curves have very similar heights, espe- 
cially for r — 2 and 3/2, away from the resonances (v x ) 
is significantly lower than Vd- The reason for this is that 
off resonance, p z increases less rapidly with increasing t, 
which means that the damping term, ap z , in Eq. ([6]) re- 
mains small and therefore causes less scattering-induced 
transport. 

Figure [TO1 reveals that, for the SL considered here, the 
most pronounced resonant enhancement of (v x ) occurs 
near r — 2 [as also shown in Fig. [5Jb)]. The two other 
resonances shown in Fig. [TO] near r = 1 and 3/2 weaken 
as r decreases because the corresponding increase of B 
reduces the distance that the electrons travel along the x 
axis in a given time interval [compare, for example, orbits 
R2 and Rl calculated for r = 2 and 1 respectively in Fig. 
He)]. _ 

To investigate how such resonances manifest them- 
selves in the electrical properties of a real device, we 
measured I(V) characteristics for the SL studied in Ref. 
[Til ], whose parameters are similar to those described in 
Section [TT1 Figure QT] shows the differential conductance 
G=dI/dV measured as function of r at V = 0.25 V and 
8 = 45°. In the experiments, we varied B and, for each 
value, measured I over a very small range of V near 0.25 
V, taking the derivative to find G. For each B value, we 
determined the corresponding r cx Fj B by assuming that 
F oc V and is uniform throughout the SL. The dominant 
feature in Fig. [TT]is the pronounced conductance peak at 
the r = 2 resonance (arrowed). Additional, far weaker, 
features can be seen for smaller r, but are not strong 
enough for us to link them definitively to particular r 
values. The overall shape of the experimental G(r) curve 
in Fig. [TT] is broadly similar to the (v x ) versus r plots 
shown in Fig. 1101 However, only the r — 2 conductance 
peak is clearly visible in the experimental data. This dif- 
ference between experiment and theory occurs because 
in the actual device, F is not constant throughout the 
device, as assumed in our present calculations, but varies 
with x due to the formation of charge domains [l4| . The 
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FIG. 11: Experimental values of the low-temperature (4.2 K) 
differential conductance, G = dl/dV, versus r = lub/uc oc 
F/B, measured for different values of B at fixed V = 0.25 
V and 6 = 45° for the SL described in Ref. fl4l ]. whose 
parameters are similar to those specified in Section [III Arrow 
marks r — 2 resonance. 



spatial variation of F means that the resonance condition 
is only satisfied in a small region of the SL, which weak- 
ens and blurs the resonances in G(r) [HI]. This effect 
is less pronounced for stronger resonances, at which the 
electric field is more uniform through the SL fl5j. Con- 
sequently, the r = 2 resonance is clearly revealed and 
correctly positioned in Fig. [TTJ 



VI. SUMMARY AND CONCLUSIONS. 

We have shown that both Hamiltonian and dissipa- 
tive chaos strongly affect the transport of miniband elec- 
trons in biased SLs with a tilted magnetic field. Chaos 
originates from a complex nonlinear interaction between 
Bloch and cyclotron oscillations associated, respectively, 
with the electric field and magnetic field component, 
Bcos9, along the SL axis. In the absence of dissipation, 
the system exhibits non-KAM chaos, characterized by 
the formation of intricate stochastic web patterns when 
u>b and uj c are commensurate. When weak dissipation 
is included a priori in the equations of motion, limit cy- 
cles replace the stochastic web filaments and stable orbits 
that they enclose. As the dissipation strength increases, 
the number of limit cycles and the volume of phase space 
occupied by chaotic trajectories both decrease. However, 
for the electron momentum relaxation rates measured for 
SLs used in recent experiments [IH, [l4| , chaos occurs for 
a wide range of F and B values, usually when 9 > 40° 
(see Figs. H ED. 

Two distinct bifurcation mechanisms drive the tran- 
sition to chaos: a period-doubling cascade, where the 
onset of chaos is gradual, and a boundary crisis, which 
switches chaos on abruptly. Both mechanisms operate for 
all tilt angles in the range < 9 < 90°. For both Hamil- 
tonian and dissipative models of the electron dynamics, 
the onset of chaos occurs near resonances when u>b and 
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u) c are commensurate. Such resonances cause an abrupt 
derealization of the electron trajectories, due either to 
stochastic web formation in a Hamiltonian picture, or to 
the creation of attractors corresponding to unbounded 
stable limit cycles when dissipation is included ab initio 
in the equations of motion. This delocalization produces 
large resonant peaks in the electron drift velocity, which, 
in turn, generate the strong resonant enhancement of the 
current flow observed in our recent experiments In 
the case of dissipative dynamics, the drift velocity in- 
creases near resonances because the most probable value 
of v x that the electron attains during its unbounded limit 
cycle is very close to the maximal possible value, de- 
termined by the SL parameters (Fig. [9]). Resonances 
between lob and lo c also give rise to multistability phe- 
nomena when two or more different attractors coexist in 
phase space. Multistability makes limit cycles, and the 
corresponding (v x ) values, evolve differently when F or B 
sweeps up or down through a resonance. We expect that 
this hysteresis will also have striking experimental sig- 
natures. In particular, our calculations suggest that ex- 
perimental 1(B) or I(V) curves will exhibit pronounced 
hysteresis, reflecting that shown in Fig. [TO] for electron 
drift velocities, especially at high field values. 

Finally, we emphasize that the driven harmonic oscil- 



lator equation ([6]), which describes the motion of mini- 
band electrons in a tilted magnetic field, is also of fun- 
damental importance in many other physical systems, 
for exam ple plas ma [17] . ultra-cold atoms in optical lat- 
tices (29l. l3Ct l3ll. the transmission of light through pho- 
tonic crystals [32|], turbulent flow and in understanding 
the patterns of quasi-crystals [HI]. Although the mo- 
tion of the dissipation-free harmonic oscillator driven by 
a plane wave is well understood, to our knowledge the ef- 
fects of dissipation have not previously been considered. 
Our study of the rich dynamics of the dissipative driven 
harmonic oscillator is relevant to a wide range of topics 
in physics, engineering, and applied mathematics. 
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